Lattice measurement of a s with a realistic charm quark 
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Abstract 

We report on an estimate of a s , renormalised in the MS scheme at the t and Z° mass scales, by means of lattice 
QCD. Our major improvement compared to previous lattice calculations is that, for the first time, no perturbative 
treatment at the charm threshold has been required since we have used statistical samples of gluon fields built by 
incorporating the vacuum polarisation effects of u/d, s and c sea quarks. Extracting a s in the Taylor scheme from the 



lattice measurement of the ghost-ghost-gluon vertex, we obtain ccf s (m^) = 0.1200(14) and a^(m%) = 0.339(13) 
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1. Introduction 



The recent announcement by ATLAS and CMS of their 
observation at 5 <x significance of a new particle with 
a mass around 125 GeV IH, interpreted as the Brout- 
Englert-Higgs (BEH) boson, makes even more crucial 
than before a satisfying control on theoretical inputs of 
analytical expression of the Higgs decay channels. In- 
deed, the era of precise Higgs physics (measurement of 
the couplings, ...) will certainly open soon: assessing 
the sensitivity of forthcoming detectors will be a key 
ingredient. There are different modes of Higgs boson 
production: however the gluon-gluon fusion is by far 
the dominant process, as shown in Fig. Q] Over the 
uncertainty Aa^ H ^ x of 20 - 25 % claimed at LHC 
( V* = 7 TeV), about 4 % come from the uncertainty 
6a s on a s (m^) J2]. A complementary approach of the 
a s measurement from the analysis of Deep Inelastic 
Scattering data, physics of jets, r decay and e + e~ — > 
hadrons [3J is its computation by numerical simulations. 
In the following section we will report on the work per- 
formed by the ETM Collaboration to measure a s from 
Nf — 2 + 1 + 1 gauge configurations fl. 
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Figure 1 : Prediction of the Standard Model BEH boson production in function 
of its mass for different production channels at LHC 



2. a s from numerical simulations 

In the past years tremendous progresses have been made 
by the lattice community to perform simulations that are 
closer to the physical point. It means including more 
and more quark species in the sea, u/d quarks (Nf = 
2), then the strange (Nf =2+1) and even the charm 
(Nf = 2 + 1 + 1) since a couple of years. Pion masses 
~ 250 MeV are now common and several collaborations 
are even able to simulate a real pion, either in a small 
volume (PACS-CS Collaboration) ^ or using a quark 
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regularisation with a rather aggressive cut-off of the UV 
regime (BMW Collaboration) [6]. Discretisation errors 
are kept under control by considering lattice spacings 
a smaller than 0. 1 fm and lattice extensions L are such 
that Lm„ > 3.5 to get rid of finite size effects. We have 
collected in Fig. [2] the simulation points performed by 
the lattice community using different quark and gluon 
regularisations. 
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Figure 2: Simulation points obtained by the lattice community in the plane (a, 



Green function [11], that we will discuss in more de- 
tails. 

The starting point is to consider the bare gluon and ghost 
propagators in Landau gauge: 



(c (2) ):y,A) 



G(p 2 ,A) 
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(F (2) ) (p 2 ,A) = -5 ab 



F(p\K) 



Choosing a MOM scheme, the renormalised dressing 
functions Gr and Fr, defined by 



Ggtf.tn = limZ 3 -V,A)G0AA), 

A— >co 



F R {p\S) = limZjV,A)F(p z ,A), 

A— >oo 

read Gr(p 2 , p 2 ) = Fr(p 2 ,p 2 ) = 1. The amputated ghost- 
gluon vertex is given by 

a q-k 

ff c (-q,k;q-k) = - O - 



F(p 2 ) 




Figure 3: Raw and cured data of the ghost propagator dressing function 

The more vacuum polarisation effects are incorporated 
in the Monte-Carlo sample, the more reliable any result 
on a s is. Several methods are proposed in the litera- 
ture to extract it: for instance analysing the static quark 
potential V(r) at short distance 17[], comparing the mo- 
ments of charmonium 2-pt correlation function with a 
perturbative formula after an extrapolation to the con- 
tinuum limit Ja], integrating the (3 function at discrete 
points in a finite volume renormalisation scheme, for in- 
stance in the Schrodinger Functional scheme [9] or fit- 
ting 3-gluons amputated Green functions in the frame- 
work of Operator Product Expansion (OPE) fioll . A last 
and particularly elegant approach consists in applying 
the OPE formulae to the ghost-ghost-gluon amputated 



= ig f bc [q v Hx{q,k) 
+(q - k) v H 2 (q, k)] . 

The renormalised vertex is Fr = Z\T; with a MOM pre- 
scription it reads 



lim Z^A) (H 1 (q,k;A) + H 2 (q,k;A))\ 2 , = 1. 

A->oo ' r 



The renormalised strong coupling constant is given by 
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In the case 



Z,0j 2 ,A) 

of a zero incoming ghost momentum k — 0, we are 
in a kinematical configuration where the non renormal- 
isation theorem by Taylor 01 211 applies: H\(q,0;A) + 
H2(q,0;A) = 1 and then Z[(yu 2 ,A) = 1. The renor- 
malised coupling in the Taylor scheme reads finally 



, 2, _ SUP 2 ) 

a T (n ) = — 

4-n 



A->oo 4n 



G(p 2 ,A)F 2 (p 2 ,A). 



The main advantage of the MOM Taylor scheme is that 
there is no need to compute any 3-pt correlation func- 
tion: it is enough to extract the dressing functions of 
gluon and ghost propagators. 

We have analysed the Nf = 2+1 + 1 ensembles produced 
by the ETM Collaboration B13I1 . with bare couplings 
P = 2.1, 1.95 and 1.9 that correspond to ag=2.i ~ 0.06 
fm, (3/3=1.95 ~ 0.08 fm and 0^=1.9 ~ 0.09 fm, respec- 
tively. Pion masses are in the rang [250-325] MeV. Lan- 
dau gauge is obtained by standard methods to minimise 
A^Afi [ 14] while the ghost propagator is computed by 
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inverting the discretised Faddeev-Popov operator. How- 
ever, as the 0(4) symmetry is broken on the lattice to the 
H(4) group, getting aj from the dressing functions G 
and F is not straightforward: 0(a 2 p 2 ) and H(4) invari- 
ants artifacts, the so-called hypercubic artifacts, have to 
be properly taken into account 1 15]: 
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Figure 6: Subtraction of aj(p 2 ) from the perturbative running 
compared to generic power corrections. 



-=() 



P [4] = 2/ Pi- We have shown in Fig. [3] that a "fish- 
bone" structure, that are those hypercubic artifacts, is 
clearly present in F but curable, as also seen on the 
plot. The remaining cut-off effects are removed by 
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Figure 4: Strong coupling constant obtained after the elimination of the domi- 
nant hypercubic artifacts 
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Figure 5: Raw data of aj(p 2 ) compared to a purely perturbative running and 
OPE with power corrections 



fitting &t(p 2 ) according to the formula '5j(a 2 p 2 ') = 
a r(p 2 ) + Cdipi o 2 p 2 + 0(a 4 ). Fig. |4] illustrates the 
benefit, in term of statistical error on 0-7-, to do the cal- 
culation in Taylor scheme, as we pointed earlier in the 



We then use the OPE formalism so that we can relate 
ar(p 2 ) to a V j A {p 2 ) flql . including power corrections: 



a T (p 2 ) = 



1 + 



d_ 



with qo fixed to 10 GeV and C a combination of a Wil- 
son coefficient and a running of the gluonic operator 
(A 2 ). Eventually is expressed at N 3 LO in function 
A _ / 507 -40AM 



of A r with = exp 



The parame- 
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ters to be fitted are thus flA^jg-, g 2 T (A 2 ), d and the ratios 
of lattice spacings 0/3/01.9 obtained by imposing that the 
various curves of aj merge onto a universal one. We 
have shown in Fig|5]that the purely perturbative running 
formula does not match with the aj{p 2 ) data, adding 
the l/p 2 fits nicely with them down to p = 3.5 GeV 
while including the l/p 6 term improves our ability to 
describe them further down to the r mass scale. One 
could expect a l/p 4 power correction but Fig|6] indi- 
cates that the fit is meaningless: still we do not exclude 
that the corresponding Wilson coefficient would mimic 
an additional l/p 2 factor. Collecting in Tab. Q] A-^g-, 
g 2 {A 2 ) and the d coefficient, with the lattice spacing 
a 1.9 = 0.08612(42) fm 10, we can run of* up to the Z° 



mass scale or at the t mass scale. For the latter we obtain 
of s (m 2 T ) = 0.337(8) and of s (m 2 T ) = 0.342(10) with our 
2 estimates of A-^g-; combining both of them and adding 

in quadrature the errors we get (rf 5 (m 2 ) = 0.339(13). 
It is in very good agreement with t decay data anal- 
ysed with dispersion relations lfl7ll . iflill . as plotted in 
FigE] Running a s up to the MS scheme b quark mass 
mb, by using the yS function and the A^_ 4 parameter we 
have measured, we can match with the Nf = 5 theory: 
a^2\m 2 ) = of^ 4 (m 2 )(l + £„ c„ (a— \m b )) n ). Then 

MS V B> MS V h ' V ^ n V MS V " / 

a second running is applied up to the Z° mass scale. 
We obtain ap(m?J = 0.1198(9) and a*f{m 2 z0 ) = 
0.1203(11) with, again, our 2 estimates of A^; com- 
bining both results and adding in quadrature the errors 
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Figure 7: Comparison of a s {m T ) measured on the lattice with estimates 
based on phenomenological analyses of r decay data flal , fl7ll . 
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we get of s (m 2 z0 ) = 0.1200(14). We have shown in FiggJ 
a comparison between lattice results H, H and JIT 
DIS data B2011 . the world average quoted by the Parti- 
cle Data Group [21] (WA '12) and a world average re- 
alised by replacing the Nf — 2 + 1 lattice results by the 
Nf = 2 + 1 + 1 one (WA '12), theoretically more reli- 
able. In [3] one can find an almost exhaustive collection 
of results. A single, very precise, lattice value domi- 
nates strongly the weighted world average of a^im^); 
removing it enlarges its uncertainty. Our estimates is 
in the same ballpark as other approaches and is using a 
complementary framework. 

Table 1 : Fit parameters of ajip 1 ) analysed by means of OPE. 

: 4 , , I T~ 1 . , . ~> . n ~T . t t " I ,1 /ft ,„ Z~Z I 



A-L (MeV) 

MS v ' 


g 2 (q 2 )(A 2 ) R 2 (GeV 2 ) 
% 


d 116 (GeV) 


a(m z ) 


316(13) 


4.5(4) 




0.1198(9) 


324(17) 


3.8(1.0) 


1.72(3) 


0.1203(11) 



3. Conclusions 

We have reported on the first measurement of a s from 
lattice simulations taking into account the vacuum po- 
larisation effects by charm quark in the so-called N r = 
2+1 + 1 theory. The main benefit of our set-up is that 
there is no perturbative treatment at the charm thresh- 
old. We have used the OPE formalism to analyse gluon 
and ghost propagators to extract the strong coupling in 
the MOM Taylor scheme. We have taken care of the 
hypercubic artifacts and included power corrections in 
the OPE, that cannot be neglected. An on-going project 
is to study whether other Green functions (3-gluon ver- 
tex, quark propagator,...) present the same feature: our 
extraction of A^g- presented here will help us to reduce 
the uncertainty on the fits of those Green functions. 
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